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We present a construction of non-equilibrium steady states in one-dimensional quantum 
critical systems carrying energy and charge fluxes. This construction is based on a scattering 
approach within a real-time hamiltonian reservoir formulation. Using conformal theory tech- 
niques, we prove convergence towards steady states at large time and we argue for their universal 
character in the low temperature regime. We compute the exact large deviation functions for 
both energy and charge transfers, which encode for the quantum and statistical fluctuations 
of these transfers at large time. They are universal, depending only on fundamental constants 
(hyks), on the central charge and on the external parameters such as the temperatures or the 
chemical potentials, and they satisfy fluctuation relations. A key point consists in relating the 
derivatives of these functions to the linear response functions but at complex shifted external 
parameters. 

February 2013 

1 Introduction 

Triggered by recent experimental advances, the study of quantum systems out of equilibrium 
is a research line of high current interest, both theoretically and experimentally. Theoretically, 
this subject is made rather difficult by the diversity of situations it a priori purports to cover, 
and by the relative lack of unifying framework extending that of thermodynamics and statistical 
physics for equilibrium systems. Situations of particular interest, which have the potential of 
being amenable to a full understanding, are those where steady currents of local quantities 
exist: steady flows of energy, charge, particles, etc. In these situations, external forces, if any, 
are constant in time, yet the system is not at equilibrium because there is a permanent creation of 
entropy. Powerful methods, based on studying fluctuations at large times and their relations with 
other physical characteristics, have been developed for classical steady states out of equilibrium, 
and it is hoped that studying fluctuation relations [12] will lead to an understanding of the 
fundamental principles of the non-equilibrium physics of steady states (and perhaps of more 
general non-equilibrium physics). 

For quantum systems, the study of steady states, and in particular of fluctuations of the large- 
time energy /charge/particle transfer, is still very much to be developed. In these situations, the 
very definition of the higher moments of large-time transfer requires the choice of a measurement 
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protocol: direct or indirect, von Neumann measurements with full wave- function collapses or 
coupling to an external measuring device with continuous partial collapse. A priori, results may 
depend on this choice. Further, as for classical systems, the construction of the steady state itself 
necessitates an appropriate set-up: external forces or coupling to external reservoirs, the latter 
being either with Caldeira-Leggett-like models [11], infinite-time hamiltonian evolution of infinite 
subsystems (hamiltonian reservoirs) [23], or non-hamiltonian evolution from integrating out the 
reservoirs (Lindblad approach) [26]. Perhaps most importantly, quantum effects themselves 
may modify our intuition on non-equilibrium physics. The full interplay between these effects 
and non-equilibrium physics is expected to be seen at quantum critical points: there, quantum 
fluctuations affect the system macroscopically, so that the dynamical transfer processes should 
be very much influenced by entanglement and interference. 

In the present paper we study non-equilibrium steady states of energy and charge transfer 
of one-dimensional systems at quantum critical points with dynamical exponent z = 1, in 
their low-temperature universal regime. Using the formulation of hamiltonian reservoirs we 
prove convergence to non-equilibrium steady states; we express these states in a precise fashion 
through a scattering-matrix formalism in the context of conformal field theory (CFT) and explain 
their universality; and we obtain exact results for all higher cumulants of large-time transfer of 
both energy and charge. The set of all higher cumulants is a very important characteristic 
of a non-equilibrium steady state, sometimes referred to as the full counting statistics (even 
though the term "counting" is more natural in the cases of countable quantities like charges 
and particles). We gather the cumulants into an exact generating function, which we will refer 
to as the large deviation function (by a slight abuse of language). From this, we verify that 
an exact fluctuation relation holds. We also show that the results are largely independent of 
the measurement protocol, and we observe that the only characteristics of the universality class 
needed is the central charge. Finally, we show that there are (at least) two fundamental time 
scales involved in the approach to the steady state: a short, microscopic (non-universal) time 
scale, and a larger one related to the speed of propagation of universal excitations. Some of 
these results for energy transfer were announced in [3], and the present paper provides the 
details of the calculations, including a description of the algebraic framework for constructing 
the non-equilibrium steady states, as well as the application of similar ideas to charge transfer. 

Many previous works have studied related quantities in different, less universal situations. 
The first exact results for the large deviation function of charge transfer in quantum systems, 
the so-called Lesovik-Levitov formula, was obtained in free fermion systems in the important 
work [25j . In other works, results have been obtained for critical free boson systems (Luttinger 
liquids) [15] , as well as in some integrable models of quantum field theory |23l [T7] . In the original 
work [25], the measurement method is indirect, the system being coupled to an external spin 
which is itself quantum mechanically measured. Other works |22[ H] use direct two-time von 
Neumann measurement protocols under various detailed prescriptions; the equivalence between 
all these protocols in the case of charge transfer in free fermion systems is a consequence of these 
works. Note that sometimes the measurement method is actually implicit, for instance with a 
physical picture of scattering particles [23] . 

Concerning energy transfer in quantum systems, the first result for large deviation functions 
was obtained, to our knowledge, in the case of a chain of harmonic oscillators in [30]. Other 
studies of non-equilibrium energy currents in the XY model gave the exact steady state and the 
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exact current [2], based on the C*-algebra mathematical development of [29], as well as some 
aspects of spin correlation functions pQ. In general massive quantum field theory, the exact 
energy-transfer steady state was described in terms of asymptotic states in [14] . Further, works 
in preparation provide various extensions of the present setting: to star-graph configuration, to 
the quantum Ising model as well as to general massive models [32] . 

Our construction is the first application of general CFT principles to the description of flows 
in steady states out of equilibrium, and the results provide the most universal critical large 
deviation formula. Since CFT describes large classes of models with interaction (in particular, 
when the central charge c is not 1/2 or 1), these are also the first results for energy current 
and large deviation function beyond models of free fermions or free bosons. With respect to the 
particular construction of the steady state used, it has been observed that it is a very nontrivial 
problem to fully describe hamiltonian reservoirs in interacting models (see for instance [6]). 
We believe our results are the first explicit interacting-model results for energy transfers in the 
context of hamiltonian reservoirs. Finally, our results clarify the universality of non-equilibrium 
steady states, and the time scales involved in its approach. 

The paper is organized as follows. After a summary of the main results in Section we 
present heuristically our set-up and the construction of the non-equilibrium steady state in 
Section [3] The formal but precise algebraic framework is detailed in Section [H This allows us 
to specify in Section [5] the non-equilibrium steady state from the construction of an appropriate 
5-matrix. Formula for the large deviation of functions for energy and charge transfer are proved 
in Sections [6] and [3 A brief discussion for our results conclude the manuscript. 

2 Overview of results 

A standard way of maintaining systems out of equilibrium consists in putting them in contact 
with two or more different reservoirs, at different temperatures or different chemical potentials 
for instance. Besides the question as to how to represent the reservoirs, natural questions that 
one may then try to answer are related to large time behaviour of such systems: Do they 
reach steady states in the long run? What are these steady states and how do the systems 
approach them? These are usually difficult questions and tools from C* algebra or quantum 
dynamical systems [1'6\ |2"S] have been developed to deal with them. Let us restrict ourselves 
to the peculiar framework - general and powerful enough at least for one-dimension systems - 
of integrable systems; these we understand as including the universal behaviours of quantum 
critical points, described by CFT. We may then wonder whether, and how, we could use tools 
from (equilibrium) integrable models to compute exact properties of non-equilibrium steady 
states. Putting an integrable system in contact with standard reservoirs (like Caldeira-Leggett 
reservoirs [11] ) usually breaks integrability, although a few exceptional cases preserving some 
integrability structures are known [20]. Considering quantum quenches [7], which consists in 
putting a system out of equilibrium by changing abruptly some external parameters, provides a 
way to circumvent this difficulty. In these cases, and in particular for local quenches in the Ising 
chain in transverse magnetic fields [8] , it has been proved that expectations of local observables 
reach steady values, and furthermore, that the system state, viewed as a dual form on the algebra 
of local observables, converges towards a generalised Gibbs ensemble. This may seem surprising 
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since the dynamics is unitary (there is no external reservoir and the system is closed). The 
key point is however that convergence only holds at infinite volume and only applies to local 
observables, so that the system is its own bath (i.e. the part of the system far away from the 
local observable locations play the role of effective baths). Note that non-equilibrium quantum 
quenches are asymptotically generalised Gibbs states, which do not support any fluxes, or energy 
or charge transfers. 

So we choose another framework, preserving integrability, which consists in preparing dif- 
ferent copies of some large quantum systems at different thermal equilibrium, say at different 
temperature and different chemical potential, and then in gluing them at a contact point. As 
soon as the contact is established, there is going to be energy or charge transfers from one part 
of the total system to another. Waiting long enough, the steady state should be established 
- this is a "real-time" construction of the steady state. Similarly to quantum quenches, and 
since transfers are localised around the contact points, the regions of the different subsystems 
far away from the contact points play the role of effective reservoirs. Reservoirs represented in 
this way using hamiltonian evolution of large subsystems at different temperatures or chemical 
potentials are usually referred to as hamiltonian reservoirs [24 ^ I29 |. 12]; they were first successfully 
used for heat transfer in |31j in the context of a classical chain of harmonic oscillators. Note 
that we may view this construction as representing what happens in the system coupled to two 
large reservoirs a finite distance apart after having zoomed over a bulk region of the quantum 
system. Under such a zoom, the reservoirs have been scaled away at infinity while preserving 
the presence of steady currents. 

Specifically, we deal with two identical one-dimensional gapless quantum systems with dy- 
namical exponent z = 1, whose low energy sectors is described by isomorphic CFTs. We prepare 
them at respective temperatures T^ r . If these CFTs have u{\) conserved charges, we also pre- 
pare them at different chemical potentials p^ r . We then connect them through a unique contact 
point, and wait long enough for the steady state to be reached. Energy and charge transfers then 
take place across the contact point, and these are steady in any finite-size observation domain. 
Gluing together two isomorphic CFTs may be viewed as modifying the defect at the contact 
point in a unique bulk CFT, from a purely reflecting defect to the absence of a defect (in fact, 
our results hold with the presence of any topological defect instead of the absence of a defect). 

First we construct the steady state, denoted cj s tat> using the algebraic tools of conformal field 
theory and vertex operator algebra. From large-time evolution with hamiltonian reservoirs, we 
show that the result can be expressed in the framework of scattering matrix, whose main input 
is to express the steady density matrix p s t a t ; dual to the steady state w s tati in a simple form, 

Pstat = S p Q S 1 , 

with S the scattering matrix and p Q the initial density matrix (note that in the following, /j sta t 
will not be defined with mathematical precision, but co s ta,t = w o Ad(5'~ 1 ) will). 

The next key point is that the steady state factorizes onto right/left movers: right (resp. left) 
movers are collectively thermalized at temperature 7} (resp. T r ), and chemical potential pi 
(resp. /j, r ) if any. That is, 

Wstat = cr AiW <g) crg ri(UT .. 

where <Tg jAt are Gibbs states on the right/left movers with f3^ r = l/(fes?] jr ), and ks is the 
Boltzmann constant. This can intuitively be understood as a consequence of the fact that 
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right /left movers have been asymptotically prepared in the far away left /right parts of the 
subsystem which effectively are at respective temperatures TJ r and chemical potentials /^ r . It 
technically follows from the fact that right/left movers have simple scattering in conformal field 
theories, even in presence of defects, and we explicitly compute the scattering matrix. 

We then use this construction to compute the large deviation function for energy and charge 
transfer. Let AtE and AjQ be respectively the energy and charge transferred during a time 
duration t in the steady regime. Then the large deviation function F(X, u) is defined by 

F(A,i;;/V,/v) : = lim 7logu; stat [ e iAA ^ e iuA ^]. 

t— >oo t 

The large deviation function codes for all cumulants of the energy and charge transferred. 

This definition is approximate, since, as we mentioned, some care is required to make precise 
the way the energy or charge transfer is quantum mechanically measured. We will consider both 
indirect, and direct two-time von Neumann measurement protocols, see Section [6j In the case 
of a direct measurement protocol, one measures the quantity to be transferred at two different 
times, and takes for AfE and AtQ the difference of the values obtained. There is an ambiguity: 
the first measurement time may be at the "contact time" , when the hamiltonian reservoirs are 
connected, or in the established steady state itself, long after the contact time. We prove that 
the result is the same for the indirect protocol, and the direct two-time measurement protocol 
with first measurement at the contact time. In fact our results also holds for the charge transfer 
statistics, but not the energy transfer statistics (contrary to the claim made in [3]), in the case 
of the two-time measurement protocol in the established steady state. 

We find a universal formula for the large deviation function, which only depend on universal 
constants (H, ks) and T\ r and \i\^ r . Namely 

F(X, u; A,r, W,r) = /(A, u; fam) + /(-A, -u; /3 r , /i r ), 
with, for unitary theories, 

J{A,v,P,») uh \ (3 _ iX p) + 8 7rh((3-i\) 8irh' 

The fact that the large deviation function F decomposes as the sum of two terms, each re- 
spectively associated to the left/right temperatures and chemical potentials, and thus to the 
right/left movers, is a direct consequence of the factorisation of the steady state. The function 
/ is a large deviation function in the chiral theory, and as such it is computable using data from 
the Virasoro or vertex operator algebras. 

The above large deviation function satisfies the fluctuation relation: 

F(i(f3 r - Pi) - \i{PlM - /3 r }l r ) - V\ /3i jT , fjL^ r ) = F(X, u; Pi jT , ni jT ). 

When chemical potentials are absent, it reduces to that which we gave in [3], in particular 
the average energy current being 



4 This is a slight abuse of language: the usual definition of the large-deviation function is as the Legendre 
transform of this large-time cumulant generating function. 

5 In the case c = 1 and T r = 0, this formula bears similarities with the Stefan-Boltzmann law for the energy 
radiated by a thermal black body, see [10]. We thank J. Cardy for pointing out this analogy. 



5 



The formula for F shows that energy and charge transfers are correlated in presence of both 
temperatures and chemical potential differences. This is not surprising if we think of these 
transfers as emanating from transfers of particles or energy and charge carrier quanta. In 
particular, a difference of chemical potentials but at identical temperatures induces non zero 
energy and charge currents, respectively: 

(Je)a/3=0A^0 = (Mi - Mr)> (Jq) A/3=0,A»¥=0 = ~ Mr)- 

Notice that in this case charge transfer fluctuations are gaussian. 

Observe also that the case where /3ini = /3 r fi r is particular from the viewpoint of the energy 
fluctuations, in that it corresponds to the case m = fi r = up to a shift of the central charge: 

p/ x n a x, _ c*tt ( iX iX \ 3 X 2 

Besides the factorisation property of the steady state, the main point in the derivation 
consists in relating the derivative of the (chiral) large deviation function / to the energy f) and 
charge j current one point functions but at shifted temperature and chemical potential. Namely, 

d 

-«^t/(A,^;^,m) = ° B jx [^(0)]) 

OA H ' P — iX 

d 

-^/(A,^;/3,M) = ^ A) »[j(0)j. 

This property is one of the consequences of PT-symetry on large deviation functions as we shall 
explain in [5]. 



3 Heuristics: non-equilibrium steady states in CFT 

The complete description of the non-equilibrium steady state of energy flow in CFT will be 
provided in Section based on the algebraic setup of CFT developed in Section HJ Here we 
give the heuristic idea behind it, starting with a general physical description of the setup and 
obtaining the steady state from standard and simple notions of CFT. 



3.1 Physical description 

We are aiming at describing non-equilibrium states in gapless (critical) one-dimensional quantum 
systems with dynamical exponent z = 1, in a setup involving hamiltonian reservoirs. The state 
is obtained by coupling two such systems which have been independently prepared at different 
equilibria; in this section we will only consider different temperatures, but the generalization to 
different chemical potentials is immediate and will be done in Sections [7] and (8) For definiteness, 
we take two identical copies of a gapless quantum system, each of length R/2, respectively 
defined on intervals [—R/2,0] and [0,R/2], and prepared at respective temperatures T\ and T r . 
We then imagine taking R large enough, and connecting the copies through the origin at a large 
negative time, say —t Q , in such a way that, after unitary evolution, the state of the coupled 
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system at time or later is stationary in any finite observation domain around the contact 
point. This steady state is out of equilibrium, and, as one can expect and as we will show, there 
is an energy flow from the high to the low temperature regions. The steady flow takes place 
on a domain of size of order Vft , where Vf is the typical velocity of elementary excitations. 
This is because, intuitively, it is excitations emitted from the contact point that allow the local 
state to change from equilibrium to non-equilibrium steady state. A similar physical picture 
tells us that for the system state to be steady, Vft Q has to be much smaller than the system 
size R. Otherwise the excitations are going to bounce back on the system boundary walls, thus 
modifying further the local states; at large times, oscillations would arise. Hence, if i is the 
observation length (the distance to the contact point up to which local observations are being 
made), the steady state is obtained in the limit R S> Vft Q 3> i. Keeping i fixed and finite, this 
means that the steady state is mathematically defined by the limits R — > oo and then t Q — > oo in 
that order H. The extreme left and right parts away from the steady-state domain, at a distance 
much greater than Vft from the contact point, serve as effective thermal reservoirs, each at its 
own temperature T\ r . Thanks to the energy flow between them, it is a simple matter to see that 
the total thermodynamic entropy of the two asymptotic regions increases. 

It is well known that conformal field theory (CFT) describes low-energy equilibrium be- 
haviours of z = 1 gapless systems (this is not proven but verified in many ways). This is because 
CFT provides their full asymptotic low-energy sector, arising from the states whose energy dif- 
ferences to the ground state are much lower than microscopic coupling energies D. CFT then 
gives a description of the low-energy thermodynamics, as well as the large-distance correlations 
of local observables in low-energy states, including, at least it is believed, dynamical correlations. 
It is important to note that these low-energy behaviours are universal, largely independent of 
the microscopic structure; CFT can be seen as giving the full quantum mechanics of universal 
degrees of freedom. It is natural to think that the non-equilibrium steady state of gapless critical 
systems, obtained from hamiltonian reservoirs as above, will be described similarly, in the region 
T\ r <C D, by a hamiltonian reservoirs setup in CFT: using a CFT state for the two identical 
copies initially equilibrated at different temperatures, and using the CFT unitary evolution from 
this state. Hence our CFT results should indeed give predictions for gapless systems. We will 
provide conditions for this to be valid in Subsection 15.21 

3.2 CFT dynamics 

Before being connected, the low-energy sector of the two copies are described by isomorphic 
CFT models with central charge c. After being connected, the total system is again described 
by a CFT model with identical central charge, but on a segment twice as large. The connection 
made at the time — t a can be interpreted in CFT as a change of the point defect at the origin: 
before contact the defect is totally reflecting, while after contact it is totally transmitting. 

Let us recall here a few standard results of CFT. In the bulk, away from boundaries and 
defects, the energy h and momentum p densities decompose as 





In much of the following we set Vf = 1, H = 1, fcs = 1 
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with h± the chiral components (right- and left-moving respectively), (dt ± d x )h± = 0. This 
implies local energy conservation, dth + d x p = 0. 

Assuming that there are no extra degrees of freedom localized at the boundaries or at the 
contact point, total, but not local, energy conservation holds with boundaries or defects. The 
boundary and defect conditions we may impose are restricted by this energy conservation. The 
total energy is the sum of the energy of the left and right parts of the system, that is 

r0~ rR/2 

H(t) := / dxh(x,t)+ / dxh(x,t). 

Demanding energy conservation, that is demanding 4iH(t) = 0, and using the bulk local energy 
conservation law dth + d x p = we get: 

= p(-R/2, t) - [p((T, i) - p(0 + ,t)] - P (R/2, t). 

The topology of our space (a single interval from —R/2 to R/2) and locality of the densities then 
forces us to impose p(—R/2) = = p(R/2) and p(0 + ,t) = p(0~ ,t). The former two conditions 
are reflecting boundary conditions at the far ends of the two subsystems, namely 

h+(±R/2,t) = h-(±R/2,t), 

which translate into the usual CFT preserving boundary conditions [9]. The latter condition 
yields the continuity of the momentum density at the contact point, that is 

h + (0~,t) - h-(0~,t) = h + (0 + ,t) - /i_(0 + ,t). (2) 

The dynamics of the uncoupled subsystems and that of the full coupled system correspond to 
two possible ways of fulfilling this condition. 

Let H Q be the hamiltonian before contact. For clarity, we will denote by h±(x,t) the time- 
evolved densities h±{x) under the i? -dynamics. The ff -dynamics corresponds to ensuring ([2]) 
by imposing reflecting boundary conditions separately on the left and right subsystems, 

h° + (0 ± ,t) = h°_(0 ± ,t). (3) 

In this case, the system splits into its two independent left and right parts; we have for instance 
h+(x,t) = h + (\x — t\) if R/2 > x,t > 0. The hamiltonian H Q is then the sum of commuting 

left and right subsystem hamiltonians, H Q = H l Q + H T with H l Q = J_d/ 2 dx h(x) and H„ = 

J^ + /2 dxh(x). 

Let H be the hamiltonian of the coupled system. We will denote by h±{x, t) the time-evolved 
densities h±(x) under the .ff-dynamics. By definition, the .ff-dynamics ensures the relation (J5]) 
by imposing continuity for the chiral hamiltonian densities, 

h±(0~,t) = h±(0 + ,t). (4) 

In this case, h±(x,t) are both continuous at the contact point, so that, for instance, h + (x,t) = 
h + (x — t) if R/2 > x, t > 0. The hamiltonian is H = f^j^/ 2 dx h(x), which formally includes not 
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only the sum of H D and H r , as did H Q , but also an extra local energy contribution at the point 
x = representing the extra link necessary in order to connect the subsystems. 

The hamiltonians H Q = H Q + H r Q and H look very similar, differing only by an extra local 
energy contribution. The correct meaning of this extra contribution is through the dynamics 
that it generates, rather than as an extra energy term (which has no clear universal meaning: the 
extra link is an energy term of order D, beyond the universal region). That is, the true difference 
between H a and H is through the defect condition imposed, eq.([3]) or eq.(jl|), corresponding to an 
abrupt change of conformally invariant defect. Before the contact time, the defect is factorizing, 
splitting the system into two parts, while after contact the defect is transmitting, letting the 
energy flow through. We will clarify this in the algebraic formalism of Section |H What we 
are describing can be thought of as a quantum quench [5], but with an initial state being the 
product of two thermal states (not just the ground state, but a statistical distribution of many 
low-energy states). 

The boundary condition (j3J) is referred to as a topological defect [28]. It ensures that the 
chiral hamiltonian densities h± are transmitted through the contact point without reflection. 
We could as well have considered connecting the two subsystems through an extra defect (like 
a quantum dot). In this case, we would need the general boundary condition ([2]), which only 
guarantees the continuity of the momentum density. It is referred to as non-topological defect. 
We shall only consider topological defect in the following. Topological defects include homoge- 
neous system with no real defect, but also certain defects making the system non-homogeneous 
but transparent to the chiral hamiltonian densities. Our results hold for topological defects in 
general. 

Note that the phase space is the same for both the H a - and i^-dynamics: it is made of the 
fields h±(x), x € M, and their descendants (discussed in Section E]). 

3.3 S-matrix and steady states 

Recall that we wish to obtain the steady state, at time 0, by preparing two identical subsystems 
at different temperatures, coupling them at time —t Q and then taking the large t a limit. The 
initial system state is the product of two Gibbs states for the left and right subsystems at 
respective temperatures T\ r . In this heuristic section we will use both the language of states as 
normalized positive linear form on the algebra of observables (to be specified later) and the dual 
language of density matrices. The connection is the usual one, which for the case of the initial 
state u , or density matrix p Q , is (where (3i >r := Tf^) 

«.[-]-*(*■■■), P.™[*-«*-™]= ^*£wy P) 

Here and in the following we use the notation n[p\ := p/Ti[p] to normalize density matrices. 

If it exists, the steady-state density matrix p s t a t, dual to the steady state w st at, is by construc- 
tion the large- R, then large-t D limit of the initial density matrix evolved with the .ff-dynamics. 
Since the initial density matrix is ff D -invariant we may write this as 

Pstat •— lim um e e p Q e e 

t — >oo R— >oo 
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This limit exists in an appropriate sense: we will show, essentially, that a similar limit defining 
w s t a t exists when acting on products of local symmetry fields of the CFT, like the local chiral 
hamiltonian densities. Note that forgetting about the defect conditions (or the single-link dif- 
ference between H a and H), the operator e~ ltoH e ltoHo is naively the identity. However, each 
exponential factor encodes the dynamics on the algebra of observables, and as such, both factors 
are distinct. By construction, it is clear that if the limit exists, then p sta t is invariant under the 
-ff-dynamics. Let the 5-matrix be defined formally by 

S := lim lim fT itoH e itoH ° . 

It is plain that the S-matrix then (formally) intertwines the initial and steady density matrices 
so that the steady state may (formally) be written as 

Pstat = S poS^ 1 . (6) 

By duality, steady correlation functions are 

w s tat[ JJ/i+^j) • JJ/t_(y fe )] := u [ YlS' 1 h + (xj) S ■ ]J S~ l h-(y k )S], 

j k j k 

with, for any local operator O, 

S~ l OS:= lim lim e~ itoH ° e itoH O e~ itoH e lt ° H ° . (7) 

t ~ >oo R— >oo 

Let us now compute the 5-matrix action on fields. By the above definition, it first acts by a 
forward time evolution with the i?-dynamics then followed by a backward time evolution with 
the _ff -dynamics. Recall that both dynamics are chiral on operators such as h± in the bulk, 
and differ only by the defect conditions at the contact point. 

Consider first the action of S on h+{x). By the .ff-dynamics this is transformed into h + (x — 
t Q ), for any R^> t Q large enough so that there is no reflection at the left end point. If x < 0, so is 
(x — t Q ) < and the field position never crosses the contact point. When applying the reversed 
-ff -dynamics one never encounters the contact point again and the field goes back unchanged to 
its original position. Hence 5 _1 h + {x) S = h + (x) for x < 0. If x > 0, then (x — t Q ) < for t S> x 
large enough and the field has crossed the contact point under the ff-dynamics. Applying next 
the backward i? D -dynamics the field will be reflected at the contact point. Since h+ is reflected 
in under the -ff -dynamics, we get S* -1 h + (x) S = h—(—x) for x > 0, so that 

S~ 1 h + (x)S = h+(x), forx<0, (8) 
S~ 1 h + (x)S = h-(-x), forx>0. 

Similarly, exchanging simultaneously left /right movers with left/right subsystems, we get 

S~ 1 h^(x)S = h-(x), forx>0, (9) 
S- 1 h.(x)S = h + (-x), for.x<0. 

This defines the action of the 5-matrix on the hamiltonian densities. Its crucial property is that 
right movers h + always end up on the left (real negative axis) , while left-movers h- end up on the 
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right (real positive axis). Using eqs, (|8|9p and the fact that the initial density matrix p Q factorizes 
on the left and right Hilbert spaces T~L l and H 7 ", it then follows that we have factorization into 
left- and right-movers in the steady states, 

WstattnM^IlM^)] = ( i j™ o w o[II /l - s (^)( _ l x jl)]) ( J^^II^-sWtel)])' 

j k j k 

with s(x) := sign(x). Identifying h±(— x) with h^(x) under uj thanks to the i? -defect condi- 
tions, we get 

c*«[nM*i) n h -(vi°)] = tn ^^)] °& m > ( io ) 

j k j k 

with (7^3 a Gibbs state at temperature on the chiral operator f) representing the hamiltonian 
densities, cf eq. (|i'2|) . This factorization and its derivation is slightly different but similar to that 
presented in [3]. 

Note that if T t = T r = T =: /3" 1 , then the density matrix dual to w sta t is e _/3i/ /Tr (e _/3H ): 
we recover the equilibrium Maxwell-Boltzmann distribution corresponding to the H dynamics 
at temperature T. We may write H = H + + H_, where H + performs the time evolution of 
right-movers (and commutes with left-movers) , and if _ that of left-movers (and commutes with 
right-movers). Then it is clear that the stationary state can be represented by 

u atat [---]=Tr(n[e-^-^ H -]---) (11) 

This is in agreement with the general understanding that in quantum field theory, the steady- 
state density matrix is an exponential of this form where H + measures the total energy of 
asymptotic right-movers, and H- that of asymptotic left-movers [HUH]. The fact that this state 
is steady, i.e. fi-invariant, is immediate to check: the i?-dynamics is left /right-moving on h±(x) 
for all x € M, and <7g is translation invariant. Tracing back the steps, the if-invariance of w s t a t 
is a consequence of the 5-matrix action (|8l9p combined with the identification h±(—e) = h^{e) 
for e — > 0, that holds under lo . Finally, note that (fit)]) holds for any chiral descendants of the 
hamiltonian densities h±(x). 

The previous argument applies to chiral hamiltonian densities and their descendants for 
any topological defect. It applies because, by definition, a topological defect is such that it 
intertwines these chiral fields on both side of the contact point. This intertwining property only 
applies to the chiral hamiltonian densities and its descendants. Commutation relations between 
the topological defect and other chiral operators may be more involved - and they are model- 
dependent details. However, if the defect is trivial, i.e. the identity so that the coupled system 
is homogeneous, the previous argument applies and the steady state measure factorizes on any 
chiral operators. In fact, it is clear that this can be generalized to any chiral current in CFT, 
for instance the U(l) current associated to charge transfer, provided the defect is topological for 
the corresponding chiral algebra. 

In the next section, we will put these heuristic arguments on a more algebraically precise 
basis. The factorization properties eq. (|10p will be instrumental. 

Finally, we note that formula ()10p and its derivation only strictly hold for local operators 
located at nonzero distances away from the contact point (i.e. all Xj ^ and all ^ 0). This is 
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first because the point x = 0, of course, needs a refinement of the initial state lu q . But second, 
and more importantly, the renormalization of fields at the point is incompatible with the S- 
matrix action; for instance, the descendant field /i_|_(0 + )/i+(0 - ) is well defined under u Q but not 
under cj s t a ti as an extra renormalization is necessary in order to take care of the zero distance 
between the two fields involved. The fundamental reason behind this is that the S'-matrix does 
not preserve the full operator algebra. We will provide a precise statement of this subtle point 
in Section [5j using the algebraic setup of the next section. This observation also relates to the 
applicability of the CFT steady state as a universal low-energy limit of quantum-chain steady 
states, a point which we discuss in Subsection 15.21 

4 CFT as an algebra of local fields and topological defects 

A precise formulation of the state w s t a t and of the S-matrix discussed above, as elements acting 
on an algebra of CFT fields, necessitates a precise formulation of certain CFT notions. Here 
we introduce such a formulation based on the description of CFT in terms of vertex operator 
algebras, combined with algebraic-QFT-like ideas of local fields. We will provide a description 
where the limit R — > oo has already been taken. 

From general QFT notions, we know that local fields cannot be multiplied if they are at 
coincident positions^. Hence we use the notion (found, e.g., in the theory of unbounded opera- 
tors) of partial algebras: the multiplication exists if and only if the positions of local fields are 
non-coincident, and we generate our basic space from this condition. On this space, we then 
define states as appropriate linear functionals. This description does not contain the informa- 
tion about the topology of the space on which the system lies: we may expect, for instance, 
the topology of co to be that of two semi-lines M+ and M_ (recall that R = oo), while that of 
^stat to be that of the line M. The topological information is provided by the operator product 
expansions (OPEs). These are equalities, in the weak sense with respect to the state, relating 
products of two local fields at different positions but near enough to each other, with infinite 
(graded) series of local fields at a single position. The OPEs provide the topology of the space, 
and give a nontrivial algebraic structure to our algebra (whose "associativity" is a consequence 
of the theory of vertex operator algebras). Given a topology there are smooth time-evolutions: 
one-parameter groups of automorphisms of the algebraic structure. 

We will consider a canonical algebra corresponding to the topology of the purely transmitting 
defect, with associated canonical state crp u p r and time evolution Ut where h± are right/left- 
movers, respectively, on the whole line. From this, an isomorphism <I> is constructed, which 
implements the purely reflecting defect by composition, see (j39|) and (|40p . Naturally, since the 
5-matrix intertwines the H Q and H dynamics, it will turn out to be exactly equal to We 
will see that the S'-matrix is an algebra automorphism for a weaker topology than that of either 

UJ or OAjtat- 

7 A non-associative "normal-ordered" product may be defined in CFT, but we will not need this here. 
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4.1 Vertex operator algebras and Virasoro algebra 

In the standard description of conformal field theory via vertex operator algebras (see for instance 
|16j ) , a model of CFT is characterized by a choice of a vertex operator algebra V (the symmetry 
algebra of the theory), and a choice of a module A4 for V (or for the tensor product V <g> V). 
Every vertex operator algebra contains the Virasoro algebra, which has basis L n : n £ Z and 
commutation relations 

rr l _ ( \t c ( m3 ~~ m ) x 

Pmi — { m ~ n )^m+n H ^ Om,n> 

and .M is naturally also a module for the Virasoro algebra. Further, a vertex operator algebra 
is naturally a module for itself, and M. contains V as a submodule (this is the identity module 
for V). The vertex operator algebra V represents the algebra of chiral symmetry currents, and 
includes, in particular, the "identity field" 1 £ V. 

The vertex operator map Y(-,x) is a map from V to formal series End(V)[[a;, x -1 ]] (series 
in the formal variable x with coefficients that are endomorphisms of V): Y(v,x) = x~ 3 Vj 
where Vj € End(V) are referred to as the modes of the element v. The sum is implicitly over all 
integers j £ Z. There is a particular element (the normalization is for later convenience) 

f> = -^L_ 2 1 (12) 

whose normalized modes satisfy the Virasoro algebra, (— 2ir)Y(fy,x) = ■ x~ J ~ 2 Lj. The formal 
series Y(v, x)w reproduces the operator product expansion of the two CFT fields v and w 
at a (formal) distance x. See [TU] for a full definition of vertex operator algebras, including 
commutativity and associativity. 

The algebra V is graded by Lq, so that V = YLjezV where j are Lo-eigenvalues and 
where every are finite dimensional. Hence we can consider infinite sums of elements lying 
in distinct graded subalgebras: Y2jez v ^ or v ^ We will refer to these as graded series. 
For instance, the formal series Y(v,x)w, for v,w £ V , x £ M x := R\ {0}, are graded series. We 
will denote by V the vector space V to which we linearly adjoin all graded series. Note that for 
v £ V and j £ Z, the endomorphism Vj acts on V. 

4.2 Algebras of local fields V and V 

Here and below, we fix V to be a vertex operator algebra as above. Consider elements of the 
form v(x) for v £ V and x £ R and let V(R) be the graded vector space spanned by these. For 
every x £ R, the map V — > V(R) : v *— > v(x) is linear and grading preserving, and l(x) = 1. 
An element v(x) will be interpreted as the field at positions x £ R associated to the vector 
v. We form the linear algebra V" spanned by all free products (juxtapositions) of elements in 
V(R) under the relations making 1 a unit. The space V is the subspace spanned by products 
with non-coincident positions under the additional relation of commutativity [£, £] ~ 0, that is, 
V = span{t>i(xi) • • • v n (x n ) : Xi ^ Xj for i ^ j, n > 1}/ ~. This makes V into a partial algebra. 
Connecting with the previous heuristic considerations, the physically relevant algebra (without 
considering the topology yet) will be V®V, with the interpretation that the first factor describes 
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right-movers and the second, left-movers, as is usual in CFT. In particular, the energy densities 
h±(x) are associated with the vertex operator algebra element f), eq. (|12|) . 

h+(x) := f)(z)®l, h-(x) := 1® f)(x). (13) 

Given a vector space W, we may make it a supported vector space by assigning to every 
C € a set Supp(C) of closed subsets of R, the supports of £, such that if A G Supp(C) then 
-B G Supp(C) for every BDi, and such that for every £, £ G and a G C, if ^4 G Supp(£) and 
-B G Supp(£) then AUB G Supp(£+£), and A G Supp(aC). A partial algebra A is consistent with 
its supported vector space structure if for every (, £ G A such that there are A G Supp(C) and 
B G Supp(£) satisfying A n -B = 0, there exists a product € A, and if for every A G Supp(£) 
and -B G Supp(£) then AU-B G Supp(££) if the product exists. The partial algebra V is consistent 
with the support structure generated by {x} G Supp(t>(x)) : v G V, x G R. 

We extend the algebra V" ^> t^ ext as follows. Let V(R) be the vector space V(R) to which 
we adjoin all graded series of the form X^'ez u 0*0 f° r u ^ ^ anc ^ £ £ R- This space 
is too big, and we consider a subspace I^ope and an associated supported subspace Fope(R) 
generated from V(R) as follows: if ( = Yljez v ( x ) an d £ = Ylkez w ^ k \y) are both in Fope(R) 
and if they have disjoint supports, then fcgZ y(u^, x — y)w^ k \y) is in V(R)ope and it has 
the support [y — d,A,y + oIa], oIa = max x£ A(\x — y\) for every A G Supp(C). In particular, the 
graded series Y(v,x)w(y) belongs to y(R)oPE for every v,w G V and x G R x ,y G R, with a 
support given by [y — \x\,y + \x\]. We note that our choice of supports guarantees, according to 
the theory of vertex operator algebras, that the sum fcgZ y(i)W, x — y)w <yk \y) indeed exists 

(using convergence in C). The partial algebra V cxt is the free commutative algebra that is 
consistent with this support structure. 

Finally, we put the nontrivial structure of operator product expansions: the algebra V is 
obtained by imposing on y ext the condition that the operator product expansion holds, 

v(x)w(y) - Y(v,x - y)w(y) °~ E 0, (14) 

for all v(x),w(y) G yoPE(R) that have disjoint supports. That is, V = V ext / °£ E . The fact 
that the equivalence relation (|14|) leads to an associative algebra product is a consequence of 
the commutativity and associativity properties of vertex operator algebra. Our construction 
V — > V cxt — > V naturally gives a map 

i v : V -> V 

vi(xi) • ■ ■ v n (x n ) (->■ vi(xi) • ■ • v n (x n ) (15) 

where on the right-hand side, we implicitly take the coset. The statement of vertex operator 
algebra commutativity and associativity is encoded into the statement that i v is an algebra 
isomorphism 

The construction (|15p induces maps, both denoted q v for simplicity, on linear functionals 
G V* and on endomorphisms^l G End(y) of V obtained by extending linearity and acting on 

8 An endomorphism K £ End(V) contains the information of a map -Ksupp between subsets of R such that if 
K{Q = £ then K S u PP (A) 6 Supp(£) for all A G Supp(C). 
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cosets in the natural way. In general, this process does not give rise to a linear functional on V 
or to an endomorphism of V. We will say that oj G V* is compatible with V if </ v (w) G V'*, and 
that K G End(V) is compatible with V if q y {K) G End(V). In these cases, we have 

q y {Co) =tDo(i v )" 1 , g v (K) = i v oio(iV. (16) 

If K\ and are both compatible with V, then their product is and we see that 

<7 V (Ki# 2 ) = q V (Ki)q V (K 2 ) (17) 

so that, in particular, g v preserves the algebra of ^-automorphisms compatible with V. 

If K is a ^-automorphism but is not compatible with V, then we may still make sense of the 
map K := q y (K) by acting on each element in each quotient. This maps quotients to a different 
collection of subsets, which gives a new algebra Vk '■ 

V K :=K(V). (18) 

Since K is a ^-automorphism, then as before, this construction induces an isomorphism 

±k ■ V -> V K 
vi(xi) • ■ ■ v n (x n ) 1 y vi{xi) ■ ■ ■ v n {x n ) (19) 

where on the right-hand side we implicitly take the associated subset in Vk (it is unique and it 
exists). The map K has the explicit expression 

K = i K oKo(i v y\ (20) 

Note that if K is compatible with V, then in particular = i v , and we recover the previous 
expression. 

When K is not compatible with V, then Vk corresponds to a different topology than that 
corresponding to V. This topology is implemented by the isomorphism defined in (|19p . 
and by the maps q^ on linear functionals and endomorphisms, defined again by extending and 
taking the quotient in the natural way according to Vk- Note that if the endomorphism L is 
compatible with V, then the endomorphism Lk '■= K0L0 K~ l is compatible with Vk', similarly 
if u> is compatible with V, then Cjk '■= & K~ l is compatible with Vk- Then, using equations 
similar to (fT6j) but for the case of Vk (in particular with i y K instead of i v ), and using (j20|) . we 
may express the maps q y K in terms of the maps q: 

q K (io K )=q(u)oK-\ q%(L K ) = K o q v (L) o K~ x . (21) 

We will use these concepts below in order to discuss defects. 

4.3 A Gibbs state and time evolution on V 

We now provide on V a family of linear functionals o~p, ft > 0, which are Gibbs states defined 
via the nonzero-temperature correlation functions of CFT. 
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The module A4 is endowed with an inner-product (•, •) that gives the usual correlation 
functions of CFT in the planar geometry (i.e. infinite- length chains in their ground states): 



c (!) Y{v-]_,x\) ■ ■ ■ Y(v n , x n ))l) 



(22) 



for Vj G V, xj G R C R 2 . Here L xi,...,x n maps the formal series (l, Y(vi,xi) • • • Y(y n , :c n ))l) to an 
analytic function of x±, . . . , x n (whose only singularities are poles) by analytically continuing from 
the region of convergence x\ > • • • > x n . In CFT, correlation functions at nonzero temperature T 
can be described via a map y i— > x from C x = C\{0} to the cylinder Im(x) G [0, /3), with /3 = 1/T; 
the map is given by x = (2-7rT) -1 logy. This map induces a map u (->• t)^ the image of which 
forms an isomorphic vertex operator algebra under a modified Y-map Y[-,y] = Y(y Lo -,y — 1) 
[27] . Hence, for ^e^itel: k = 1, . . . , n, we define erg G V™ by 



Cr/3[fl(xi) • • •■UnOCn)] 



., yn (l,Y( yi L °vf, yi ) 



Y(ytv^\y n )l 



(23) 



We note in particular that (see f)12|) and note the normalization of fj) 



(/?) = (?HL 



fi> VJ ' 48vr' 

By the theory of vertex operator algebras, dp is compatible with V and 

op := q v {dp) G V*. 



(24) 



(25) 



Correlation functions at nonzero temperature are also traces on appropriate modules. By 
standard arguments, CFT for a chain of a finite scaled length R/2 has a module M for a single 
copy of V describing states on a circle of circumference R. The state of nonzero temperature 
T is geometrically a torus with cycles of lengths R and (3 = 1/T, obtained by tracing over the 
module M: 



8/1, 



,yu^M (n 



e 



y k =e—z-} 



(26) 



Here the element Lo of the Virasoro algebra is involved. The infinite- R limit then gives the state 
of nonzero temperature on the line: 



lim 

R— >oo 



v yi,—,y-> 



e 



Lo 



r Hnl£fe-| • (27) 



We mention that from the data of the one-point functions op \y{xy\ (which are in fact inde- 
pendent of x), it is possible to deduce the correlation functions ap[vi(xi) ■ ■ ■v n (x n )'j by seeing 
them as giving values on the real subspace R n of functions on C n with appropriate analytic 
properties, including periodicity along the imaginary directions with period f3 and poles at col- 
liding positions of a fixed form determined by the vertex operator algebra. These are 
the conformal Ward identities. 
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There is a natural time evolution associated to the Gibbs states defined above, which pre- 
serves the state. Let : t € R be a one-parameter group of ^/-automorphisms defined by 

U?( v (x)) = v(x -t), ueVjxER. (28) 
These automorphisms are compatible with V and we define 

■= q y (U?) € End(V) (29) 

which, thanks to (|17p . form a one-parameter group of automorphisms of V. Physically, the 
compatibility indicates that represents a time evolution in agreement with the topology on 
the line, here for a single copy of V, and the definition (j28[) means that this single copy is 
right-moving. These automorphisms preserve the Gibbs state or defined above: or o Uf = a p. 

4.4 Canonical topology, states and time evolution on the physical algebra 

v = v®v 

The physically relevant algebra is V <g> V. This is naturally a supported vector space with a 
slightly extended notion of support, where supports are generated in the natural way from those 
of V and lie on the disjoint union K Udi s j K- We now consider on V (8> V a canonical topology, 
canonical states and a canonical time evolution. 

The canonical topology corresponds to performing the construction (|15j) on both factors, 
leading to maps q = q v ® q v and i = i v <8> i v . We will denote the corresponding algebra by 

y := V®V. 

Canonical states are defined by factorization: 

a /3i,/3 2 = <rpi ® o-fc = q{ap u/32 ), fix, fa > (30) 

where r 2 := dp x ® cr/3 2 (these are compatible with V). Finally, the canonical time evolution 
is the one-parameter group Ut ■ t G M of V (8> ^-automorphisms defined by 

U t = C/ t v Ul t = q(p t ) (31) 

where Ut '■= U^ ® U^ t (compatible with V). In terms of our heuristic discussion above, where 
the hamiltonian H is involved, the time evolution Ut(---) is interpreted as the conjugation 
e lHt ■ ■ -e~ lHt ; it physically represents the evolution on the line, where in V = V (8> V (hence 
also in its restriction V <3 V), the first factor is right-moving and the second is left-moving. 
In particular, this time evolution, along with the correspondence (|13p . implies the usual time 
evolution for right- and left-moving energy densities, 

h±(x,t):=U t (h±(x)) = h±(xTt). (32) 

The automorphisms (|3T|) preserve the linear functional (yp u p 2 , 

<T/3 1 ,fc°Ut = (Tp u p 2 . (33) 
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The topology represented by V puts both right- and left-movers on separate lines, so that 
the full quantum system, containing both types of movers, is on the line. Note however that 
0/3 1|( a 2 is not in general a physical Gibbs state for the corresponding quantum system, except 
when pi = fa = 13, where we then have a state of nonzero temperature f3~ l on the line (which 
we do not make use of in this paper) . 

4.5 Boundaries and topological defects 

In our heuristic discussion we introduced two types of defects: a reflecting defect at x = 
and a topological defect, effectively the absence of any defect from the viewpoint of symmetry 
currents. These two types of defects were associated with particular time evolutions. Also, on 
systems with the former type of defect we put a product of Gibbs states where both halves are 
thermalized independently at temperatures (3i and f3 r . On systems with the latter type of defect, 
we had instead a non-equilibrium steady state. We now describe how to implement these defect 
conditions on the local field algebra V ® V, both from the viewpoint of the correct topology and 
of the time evolution, and we explain how to implement the initial product of Gibbs state on 
the former type of defect. We will describe the non-equilibrium steady state in the next section. 

The main idea is to intertwine the physical degrees of freedom in the topology-less algebra 
V <S> V in such a way that the topology q, the canonical time evolution Ut and the canonical 
state (Tfa^fc §i ve r ^ se to the physical topology, time evolution and state. The intertwining of the 
physical degrees of freedom will be implemented by a V <8> ^-automorphism <3?, 

§ :V ®V ®V. 



As in psp. there is a natural extension 

$ := g($) : V -»• V$ (34) 

which provides the algebra V$, with associated maps q$ and isomorphism i$ : V ® V — » V$ (see 
(|19p ). We have, explicitly, <3? = i$o$oi _1 (see (|20p ). This means that the map <E> -1 intertwines 
both the local degrees of freedom and the topological information of the canonical algebra V in 
order to give the physical algebra V$ 

Given <J>, the physical time evolution and physical state are those obtained by extending 
(C/t)$ := $o[/ t o and (5' ( a 1)j a 2 )$ := ^/3 1 ,/3 2 ° < £ -1 > the latter for appropriate $\ and /?2- 
Following (|2ip . these extensions are respectively 

g*((l/t)*) = $ ° U t o and g*((cr i 9 1 , j 9 a )*) = a^fo ° (35) 

In general, $ is not compatible with the canonical topology, so that V$ represents a different 
topology than that on the (double) line MUdisjlR- If it is, then $ = io^oi" 1 is a V-isomorphism, 
so that V$ = V. In this case, we recover the topology on the line. 

For describing the system with a topological defect, we simply choose $ = id, in which case 
$ is a (trivial) automorphism. 

In order to describe the reflecting defect, we set $ = defined by 

v(x) <g> 1 v(x) ® 1 (x < 0) 

~ _ h-> l®u(-x) (x > 0) , . 

°" l®v(x) i — ^ u(-x) (2) 1 (x < 0) ^ ' 

1® -y(x) i — y l®v(x) {x > 0) 
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for every v 6 V, with supports transforming as 

: (A. U A+) U disj (5_ U 5+) m- (A_ U (-£_)) U disj ((-A+) U B+) (37) 

where .A_ C (— oc, 0], A + C (0, oo), C (— oo, 0) and B + C [0, oo). This has, physically, 
the effect of "unfolding" the left-hand side into right-moving currents, and the right-hand sides 
into left-moving currents (or, of folding the right-moving line into the left-hand side, and the 
left-moving line into the right-hand side). Writing V <S> V = (V<o <8> V>o) <8> (V<o <8> V>o) in an 
obvious notation, the V <g> ^/-automorphism $ Q is essentially a permutation of the middle two 
tensorial factors (with appropriate permutation of the associated supports). Note that there are 
other consistent choices for the inclusion of the boundary point x = 0, and that this is largely 
irrelevant for our discussion. 

The full local field algebra in the case of a reflecting defect is then V '■= V$ D and the 
corresponding physical intertwiner 

3> := q{$ ) : V -> V Q . (38) 

We note that $ is not compatible with the canonical topology. This is physically because V ®V 
contains the OPEs representing the topology of the line, in particular there is a nontrivial OPE 
expanding h + (— x)h + {x) for x > into a graded series of fields at x; while with a cut, we do 
not expect h+(— x) and /i+(x) to interact in any way. 

The time evolution on V compatible with the reflecting defect is then 

U° := $ o U t o fcj 1 . (39) 

Thanks to the unfolding performed by according to this time evolution there is reflection 
on the boundary at x = both for fields coming from the left and fields coming from the right. 
Although U° is not a time evolution where v(x)®l is always right-moving and l<g>v(x) is always 
left-moving, still for every i / there is a t' > such that for every v,w € V and |i| < t', 
we have U°(v(x) (8 w(x)) = Ut(v(x) <S> w(x)) = v(x — t) <g) w(x + t). This can be generalized 
to other elements as long as they are supported away from the point x = 0, so that U° is 
still, locally and away from the defect, a time evolution where the first factor is right-moving 
and the second left-moving. The only difference with Ut occurs at the point x = 0. In terms 
of our heuristic discussion, this time evolution is obtained from the hamiltonian H Q , so that 

[/?(...) O e iH t . . . e -W t 

Finally, the unfolding also guarantees that the state u where left and right subsystems have 
been separately thermalized can be represented by 

u = <r A A o (40) 

By construction, this state is preserved by the time evolution U°, oj q o U° = uj . It is worth 
noting that uj factorizes on the left- and right-subalgebras, thanks to the exchange of the middle 
tensor factors produced by $ G as described above, and to the factorization of ap^fc on the first 
and second tensor factors of V (8) V. 
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5 /S-matrix and the non-equilibrium steady state 



5.1 The CFT non-equilibrium steady state 

As explained in the heuristic discussion, the real-time (hamiltonian-reservoir) construction of 
a non-equilibrium steady state with thermal flow through a topological defect is obtained by 
starting with two systems with boundaries, thermalized at different temperatures, which are 
then connected by a topological defect and evolved for an infinite time according to the time 
evolution of a system with a topological defect. The resulting steady state is a state with energy 
flow on the line, and is stationary with respect to the evolution on the line Ut- 

In the present subsection, we show that this large-time limit exists weakly on V, we establish 
the form of the resulting stationary state w s tat> stationary with respect to Ut, and we show 
that it can be obtained from the initial state oj by composing with a "scattering matrix" (or 
scattering map) S : V — > V which is an algebra isomorphism (this corresponds, in our heuristic 
discussion, to Ad(S' _1 ), whose action is given by ([7])). In particular, the scattering matrix is not 
an automorphism of the initial algebra V , which indicates that it implements a change of the 
locality structure of the algebra of observables. It is an automorphism of an algebra formed by a 
common subspace of V and V Q , denoted Vb<g)Vb, where essentially the point x = has been taken 
out altogether. But this algebra has a weaker locality structure which is not preserved by the 
time evolution Ut- The fact that the space topology preserved by the scattering matrix is weaker 
than both those of the defects is an indication that certain QFT renormalization processes are 
modified in the steady-state limit. This reflects the fact that certain high-energy states may 
occur after the steady-state limit which are not described by CFT, and that the associated 
observables do not have a steady-state limit in the CFT context. We will discuss this aspect in 
more detail in the next subsection. 

In order to describe the algebra Vq <8> Vq and in fact to discuss the large time limit generating 
the steady state, we now make the following construction. For every s£l consider the subalge- 
bra V <s of V formed by all elements supported in (i.e. possessing supports in) (-co, s), and the 
subalgebra V >s of V formed by all elements supported in (s, oo). The algebra V s := V <s ® V >s 
represents a topology where a single point has been taken out at x = s, thus separating the left- 
from the right-hand side of s; OPEs are implemented just for local fields all standing on the left 
or all standing on the right of the point x = s. Note that this is not like a reflecting defect, 
because we have only one copy of the local algebra, so there is not reflection phenomenon. The 
state erg defined in (|25p gives rise to a well-defined linear functional on V s for every s. However, 
the time-evolution operators C/ t v as defined in (|29jl are not automorphisms of V s . This is because 
right-movers will eventually go through the point x = s. Yet t/ t v provide isomorphisms between 
the algebras V s , and we have 

C/ f v : V s -> V s - U a,t€R. (41) 

Now form the space V s <8> V- s . This is a sub-algebraic structure for V. Hence the state o" i a lji g 2 
is a linear functional on it. However, the physical state representing the reflecting defect uj q in 
which we are interested is generically not. In fact, it is a linear functional on V s <8> V- s if and 
only if s = 0. This is because for s = 0, there is no OPE implemented between fields on the left 
and on the right of the defect, so oj q is well defined on the algebra. 
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We find that Vq ® Vb is a sub-algebraic structure for both V and V , and that $ as defined 
in (|58"]) is a Vo Vo-automorphism: 

$o : Vb <8> Vb ->• % ® Vb- (42) 

We note also that the limit lim^oo V s <E) V- s exists in a local sense and gives V (in particular, 
for every Q € V there exists a s' such that £ € <S> VL S for all s > s'). 
We can now express our main proposition. 

Proposition 5.1 The limit 

lim U o [U t (0] =f«rtat[C] ( 43 ) 
t— >oo 

exists and is finite for every ( € V, and gives rise to a stoto Wgtat on V given by 

Wstat=CT AA . (44) 

Proof. First, we observe that thanks to (|4"Tj) and ([3"T]) . the state ui oUt is well defined on Vt®VLf 
Hence by the discussion above, if the limit exists, it provides a state on V. Further, since uj q is 
invariant under U°, then we may evaluate 

limu; [((U t °r l oU t ) (C)]. 

It is sufficient to consider £ = i(£) for every £ 6 V" <g) V", and it is sufficient to take generating 
elements in V®V. We now evaluate the action of (U°)~ l oUt on generating elements t>(x)<g)l and 
l(2)v(x), for v G V" and x£l. Thanks to (|3ip . for every x, there exists a i large enough such that 
Ut(v(x)<S>l) is supported on (oo, 0) and Ut(l<Siv(x)) is supported on (0, oo). Then, using ([36]) and 
(|38p . we see that on these vectors $ = id. Hence, using (|39p we find that for every x there exists 
a t large enoug h such that (({7 t ° )" 1 o C^) (u(x) ® 1) = ($ o J7 t _1 o I7 t ) (u(x) ® 1) = $ (v(x) <g) 1) 
and similarly ((U°) o [M (l<gw (x)) = $ (l(g)v(a;)). Since this holds for all generating elements, 
this shows that 

lim u o [U t (0] =Kod> ) [(). 

Using (00]), this completes the proof. I 
This can be re-expressed in terms of the scattering matrix <S. 

Proposition 5.2 We have that 

S : = lim ((U?)- 1 o U t ) = $o : V -)■ V G (45) 

in i/ie sense that for every ( £ V, £/iere is a t large enough such that the left-hand side of the 
following is well defined and 

((U t °r 1 oU t )(0 = M0- (46) 

Then, 

Wstat =UJ Q oS. (47) 

Proof. This follows immediately from the discussion above and the proof of Proposition 15.11 I 
We finally note that thanks to (j4"2]l and the discussion around it, we have that S is a Vq^Vq- 
automorphism: 

S : Vo ® % -> Vq <8> Vb. (48) 
That is 5 preserves a weaker locality structure than that of the full algebras V or V . 
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5.2 Universality and scaling limit 

It is common belief, without any reasonable doubt but in most cases without formal proof, that 
1+1-dimensional CFT describes the low-energy sector of quantum chains (one-dimensional chain 
of quantum degrees of freedom connected by few-neighbour interactions) at critical points with 
dynamical exponent z = 1. That is, in such critical systems and in the thermodynamic limit, 
whenever only states of energies much smaller than the typical microscopic interaction energies 
(or, in conduction electron systems, the bandwidth) D are excited, then the resulting physics 
is correctly described by CFT. This includes the equilibrium thermodynamics as well as the 
large-distance coordinate dependence of equilibrium correlation functions of local observables. 
The description provided by a given CFT model is common to every microscopic model in a 
universality class: the low-energy sector is universal. We would like to claim that this holds as 
well out of equilibrium: under natural assumptions, the CFT description of the non-equilibrium 
steady state we have given above correctly describes non-equilibrium steady states of critical 
quantum chains (as obtained in a real-time construction with hamiltonian reservoirs) at low 
temperatures I) )T ., and it is universal. 

Namely, given an observable in the quantum chain, we make three assumptions. The first 
assumption is that the stationary limit exists for every temperatures 7} and T r . This, thanks to 
the Gell-Mann-Low theorem, implies that when 2} = T r , the limit gives the equilibrium Gibbs 
state of the connected system. The second assumption is that the low-energy sector is described 
by CFT, in agreement with common belief. The third assumption is that with 7} = T r , the 
low-temperature limit of the equilibrium average of the observable, both before contact and 
after the stationary limit has been reached, converges to the CFT average. From this, we argue 
that our CFT computations describe correctly the low-temperatures limit, T\ T <C D , of the non- 
equilibrium average. Note that this is nontrivial because the quench that we perform, which 
amounts to connecting two halves of the system, involves adding a finite number of microscopic 
links, each of energy of the order D. Note also that the low-temperature limit necessitates in 
general a renormalization of the observable - we will understand the observable as flowing with 
the scale of the temperatures in an appropriate way, and discuss the subtleties related to this. 

To get a hint on the usual argument for the convergence toward CFT at low energy, let 
us write the equilibrium quantum chain Gibbs state at temperature T = (3~ l as w* a : = 
J di/g ha (E) e"^ E w*g, where wjjj^ is the normalized state that traces over the subspace of energy 

E, and du^(E) is the energy density in the thermodynamic limit. Let O c ha be an observable on 
the chain, whose average we wish to evaluate. Provided that the energy density does not grow 
too fast at high energy E 3> D, and likewise that the trace w^fe [C c ha] does not grow too fast, 
all contributions from energies E 3> D are exponentially suppressed by the factor e~@ D in the 
expectation uf a [O c h a ] • Then, with an appropriate choice of O c ha 5 renormalized appropriately 
with the scale (3D, the usual assumption is that Wo ha [O c ha] converges towards the CFT Gibbs 
state oj as (3D — > oo, such O c h a is mapped to an appropriate CFT fields (or product thereof) 
O. 

The hamiltonian-reservoir real-time description of the non-equilibrium steady state we pre- 
sented above may be applied to the quantum chain. The initial state w* a represents a direct 
product of Gibbs states at temperatures T\ and T r for two separate half-infinite critical quantum 
chains, with hamiltonian H^ ha = J7 ; cha + if£ ha , and the evolution hamiltonian is i7 cha , where 
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the interaction links between left and right half-chains have been added in order to make a 
homogeneous infinite chain. We may represent w* a as above, but by separating the energies of 
each half 

.cha 



0J„ 



J dvt*{E h E r ) e-^~^u^ Er , (49) 

where ^^ E[ E corresponds to a trace over the Hf 1 ^, i^: ha -eigensubspace with fixed Ei,E r . 
Then, the quantum chain steady state cj cha may be written as 

[©cha] = wf* [^OhaScha] (50) 

for an appropriate family of operators O c h a ( e -g- operators supported on a finite number of sites), 
where 5 c h a — limbec e~ liHA ™ e itH ° * is the lattice S- matrix. We assume that the stationary limit 
t — )■ oo exists for all fa,/3 r . This implies that the stationary limit exists in w^j 1 ^ E [S' c ^! 1 c ha5'cha] 
independently for every Ei,E r . From Gell-Mann-Low's theorem, the 5-matrix intertwines the 
hamiltonians H Q and H, 

H c a 5* c ha = S c h a Hg a . (51) 

This can be understood from the existence of the stationary limit, formally by setting to zero 
the time-derivative: lim^oo d ^ e _4 *^ cha e*'^° ha ^ /dt = 0. This implies that S^ha maps each 

Hf ha , ff^ ha -eigenstate of eigenvalue Ei,E r to an .ff-eigenstate of eigenvalue E\ + E r . There- 
fore, for every Ei, E r , there is a subspace of the if-eigenspace of energy E = E[ + E r , such that 
the associated density is 

dv cha (E l ,E r ) = d^ ha (E l ,E r ), 
and such that the associated trace is 

<E r = ^ uEr oAd(S^). 

Hence, 

J dv cha (E h E r ) e ~^ E '-^u;fX- ( 52 ) 

Note that if du cha (E) and uj c e & represent the density and trace state for the -ff-eigenspace of 
energy E, then we have f El+Er=E dv cha {E h E^ojf^ = dv cha {E)ujf a . With fa = fa = fa we 
see that Wgtat i s the Gibbs state for the hamiltonian H at inverse temperature fa 

We now consider the low-temperature limit. The assumption that the low-temperature limit 
exists and is described by CFT at fa = fa both before the contact (i.e. with wj; ha ) and in the 
stationary limit (with w^.), and the assumption that the low-energy sector is described by 
CFT, imply that in the integrals (|49p and (|52|) . the low-temperature limit is obtained from 
the low-energy sector E\ + E r <C D. This means that at large energies, the averages under 
oo^ E and oo^ E and the energy density do not overcome the exponential suppression, and that 
the renormalization is that necessary for the limit to be finite. Hence this also holds with 
fa / (3 r : as both temperatures are lowered, only the low-energy sector is involved, so the result 
is described by CFT. It is clear that our CFT construction for u describes the low-temperatures 
limit of uj ha . However, the meaning of E[ and E r in the CFT limit of the integral (|52p is not 
a priori evident. In order to argue that the low-temperatures limit of uj^t is w stat> we need 



cha 
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to make sure that the limits of large time and of low temperatures can be interchanged. We 
argue from the idea that in the limit (|5U|) , the evolution can be seen as an adiabatic change from 
eigenstates of H Q to eigenstates of H (there is no adiabatic switching of interaction, contrary to 
the setup of the usual Gell-Mann-Low theorem, but the presence of the observable O c ha plays 
this role). Since both at the beginning and at the end of the process, only low-energy states are 
necessary to describe the average when temperatures are low, then this suggests that throughout 
the process the low-energy sector is sufficient to describe the low-temperature limit. Hence, the 
whole process can be described by CFT, and we recover our CFT construction of w stat . 

The arguments above make it clear that w s t a t describes the low-temperatures limit of the 
quantum chain steady state only for observables whose renormalization is the same before the 
contact and in the stationary state; otherwise we cannot guarantee that the process from the 
former state to the latter state corresponds to an adiabatic change, hence we cannot interchange 
the low-temperatures and large-time limit. Due to the local quench, local operators near the 
origin will have different renormalizations. For instance, the product of degrees of freedom on the 
left-hand side by degrees of freedom on the right-hand side, both near the origin, renormalizes 
differently before contact than in the stationary state. This means that we have to restrict the set 
of quantum chain observables for which our results describe the correct low-temperatures limit. 
This restriction is the quantum-chain underlying phenomenon corresponding to the reduction 
of the algebra of CFT observables from V to Vq <S> Vq . 

5.3 Approach to the steady state 

To decipher how the system state converges towards its steady states we have to look at the 
large but finite time t behavior of the correlation functions 



where the time evolution is defined by the if-dynamics. As soon as t Q is larger than the maximum 
of the distances of the field locations from the contact point, that is t Q > max[\xj\, \yk\], all the 
chiral operators h± have moved either to the left or the right subsystem since then xj — t Q < 
for all j and yt + t a > for all k. The above correlation functions then factorize on left and 
right sectors as 



and these two factors are t independent due to the iT -invariance of the uj q measure (recall 
that translating h± in the w Q -correlation functions using the ff -dynamics may require using the 
reflecting boundary conditions h±(x) = h±(—x)). 

Hence, on chiral hamiltonian densities the approach to the steady state is abrupt. It is 
instantaneous as soon as all operators have moved to their asymptotic subsystems, and this 
takes a time of the order of the maximum size of the field location domain. This is valid in the 
continuous limit. In the quantum chain model with a microscopic energy scale D, we expect two 
characteristic time scales, the large infrared time scale of the order of the size of the observation 
domain associated to an abrupt change as we just described, and a short non-universal time 
scale of order 1/D associated to exponential convergences, as seen numerically |21j . 






j 
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6 Energy flow and its large deviation function 



We now evaluate explicit physical quantities form the stationary state w s t a t : the energy current, 
and the large-time statistics of the fluctuations of energy transfer. 

6.1 The energy current 

By standard field theory arguments, the local operator whose average gives the energy current 
is the momentum density operator ([I]), 

Je(x) = p(x) = h + (x) — (53) 

with h+(x) and h-(x) the right-moving and left-moving chiral energy densities (|13|) . Thanks 
to Proposition 15.14 the stationary average is w s tat [</e( x )] = ^Pi,p r \Je{%)\ ; an d thanks to ([32]) 
and (|33|) . it is independent of x, so we can choose x = 0. Then, thanks to the factorization 
(|30|) . and using the expression (|23p . (|25|) as well as translation invariance on the plane (i.e. the 
L_i-derivative property of vertex operator algebras and the fact that L_il = 0), we have 

^stat[J E (o)] = (i,Y(f^),o)i) - (i,y(f,(*>,o)i). (54) 

Finally, using ([24]) and the fact that (1, Y((j,0)l) = 0, we find 

oo stat [J E (0)} =^(^-2-^-2) (55) 
in agreement with the expression for ( Je)a/3^o,a^=o quoted in Section [2] 

6.2 Large-time energy transfer statistics 

6.2.1 Definition, measurement protocols, proposition 

Let us refer for a moment to the heuristic discussion of Section [3] In order to evaluate the large- 
time statistics of the fluctuations of energy transfer, we need to define an operator whose time 
derivative (with respect to the -ff-evolution) is the energy current, and evaluate the statistics 
of the change this operator undergoes after a large time t. A natural candidate is the operator 
that measures half the total energy difference between the left- and right-hand side: 

After the infinite-length limit R — >• oo has been taken, the operator E does not make sense 
anymore, as it then measures the difference between infinite quantities; technically, it also falls 
outside our formalism, as it is not a local operator, and in particular it does not have a stationary 
limit. Yet, after a finite time t, its change is 

AE(t) := E(t) — E = ( J E (0,t)= [ dx (h+(-x) - h-(x)) (57) 

which is expressed in terms of local fields and which possesses an infinite-length and stationary 
limit. This suggests that the statistics of the change it undergoes after a finite time t can 
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be evaluated in the non-equilibrium stationary state using our local-field formalism. We now 
explain how this is done. 

In order for the initial part of the derivation to make sense, let us first consider a finite- 
length quantum chain (hence in particular before taking the scaling limit, in order to avoid UV 
divergences), on which we have an operator E defined as in (|56p . One has to be careful as to 
how to define the energy transfer during a time t. An experimentally relevant and standard 
way is to use an indirect measurement protocol: essentially, the current Je is linearly coupled to 
another quantum system, a "dial" that increases proportionally with time with a proportionality 
constant given by Je- Quantum measurements are performed on this dial, and from this, the 
statistics of the energy transfer is inferred. In particular, one infers the probability Pt(Ae) that 
the energy transferred after a time t be AE(t) = Ae, and defines 

Z t {\) :=^e* AAe P t (Ae). 

Ae 

One then constructs the large deviation function 

F ind (A) := limt-Mog-Z^A) (58) 

t— >oo 

This indirect-measurement setup was considered in |25| in the case of charge transfers, where 
more explanations can be found. Adapting the expression found in |25j to the case of energy 
transfer, we then expect 

F ind (A) = limf 1 log (Tr (p stat e^e^)) . (59) 

In fact a formula very similar to (|59|) can be obtained if we consider standard von Neuman 
measurements of E, instead of indirect measurements, as follows. Assume two measurements 
of E are performed: First E is measured exactly at the original contact time —t . The output 
is ei with probability Tr(P ei p), where P ei is the projector on the corresponding eigenspace. 
Then, after a long time, at time 0, E is again measured. The output is &2 with probabil- 
ity Pt (e2,e\) = Tv(P e2 e~ ltoH P ei p P ei e lotH P e2 ). Since p Q commutes with P ei , using the main 
property of projectors we can simplify this to Tr(P e2 e~ ttoH P ei p e itoH ), which we transform into 
Tr(P e2 (t )P eiPo ). Using £ e /(e)P e = f(E) and P t (Ae) = £ ei , e2:e2 _ ei=Ae Ptfa, ei), we obtain 

F vN (A) = Mr 1 log (Tr ( Po e^We"^)) . (60) 

We remark that although the operator averaged in eqs. ([59]) and (|60p appears non-local, 
the Baker-Campbell-Hausdorff (BCH) formula guarantees that only AE(t) = E(t) — E and its 
commutators with itself and with E are involved. Since the time-derivative of E(t) only involves 
an operator supported on a finite number of sites near the origin, then AE(t) is essentially 
localized near the origin, on a region of length proportional to t; whence also all commutators 
with itself and with E are. This suggests that the infinite- length and, as is needed in (|59p . sta- 
tionary limit exists; and, in accordance with our arguments of Subsection l5.2| that in the scaling 
(low-temperatures) limit it is given by evaluating the corresponding average in the stationary 
state w s t a t in CFT. In the CFT, by the BCH formula, the operator averaged in (|59p and (|60p is 
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explicitly a series of integrals of local fields supported on a region near to [— t, i]; whence indeed 
its stationary limit exists in the CFT construction itself. It is this integral expression, written 
in (|74p . that we will use as our starting point for our actual proof of Proposition 16.11 below. 
We will show the following. 

Proposition 6.1 The energy cumulant generating function f7^| ) in the non- equilibrium steady 
state w s tat of CFT is 

F(A):= FvN (A) = ^(A) = ^(^^-^^) (61) 
where /3i >r = (k B T^ T )~ l . 

Here we have put back Planck's and Boltzmann's constant. Note that the result is very universal, 
depending only on the CFT central charge and universal constants. The order 0(iX) gives the 
current w stat [ Jg (0)] = (Je)a / 9^o,a^=o quoted in Section [2] and derived in ([55]) . The function 
(f6"T|) satisfies the fluctuation relation [121 US] , 

F(i((3 r -/3 l )-X) = F(X). (62) 

That the energy transport fluctuations satisfy the fluctuation relation has been checked in the 
quantum harmonic oscillator chain [30]. As usual, the fluctuation relation relates the probabili- 
ties Pt(0) and Pt{—6) of opposite energy transfers E(t) — E = ±t6 across the interface: 

We mention here that it is possible to define a two-step measurement protocol as described 
above, but where the first measurement time occurs when the steady state has been estab- 
lished (time 0), instead of at the contact time (—t ). Similar arguments as those above lead to 
Pt(^2t e i) — Tr(P e2 e~ ttH P ei p s t a tP ei e ltH P e2 ). For finite lengths, E has a discrete spectrum and 
we may use |15] the formula J due lu ( E ~ ei ^ oc P ei (with an appropriate integration range). This 
yields, naively, an integral representation 

1 [ T 

Z t (X) o: lim — / duZ t (X,u) (63) 
r->oo r J_ T 

with 

Z t (X, u) := Tr ( p stat e<^) E e iXE ® e<-z +U ) E ) . (64) 



Here r — > oo is the inverse of the typical (vanishing) energy difference between nearest states in 
the spectrum. Many subtleties occur in a proper argument leading to (|64|) ; and further, from 
(|64p . one finds that the associated large-deviation function is not ()6ip . See appendix lAl for 
calculations pertaining to this. 
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6.2.2 Proof of Proposition 



We first explain the lines of the proof by referring to the heuristic discussion of Section [3l 
Consider the expression (fTT|) for w s t a t [O] and take 



O 



e i\{E+AE(t)) e -iXE 



(65) 



Thanks to (|57p . the difference AE(t) contains right-movers only on the left-hand side, and left- 
movers only on the right-hand side. By (|56p and the fact that both H. T Q and H act like H^ + H- 
away from the origin and on the right-hand, respectively left-hand side, this means that all 
commutators involving E and AE(t) can be obtained by replacing E i— > %{H- — H + ). Hence, 
thanks to the BCH formula, we have the factorization 



,iXE{t)-i\E] 



where 



Z t (X;f3) :=Tr n 



-ft #4 



Z t (A;#)Z t (-A;/3 r 



,%H+ -%H++iXf*dxh+(x) 



(66) 



(67) 



The argument above also means that we may do the replacement H T i— > H + and H l Q \— > H- in 
the density matrix p Q itself in ([60]) . which corresponds to p Q i-> p stat and shows the first equation 
of d6U). 
Define 

/(A;/3) := lim r 1 log Z t (X; /3). 

t— >oo 

Note that the argument of the limit in t is to be understood as a series expansion in A (a 
generating function). Taking the A-derivative on this series, we have 







d(iX) 



/(A;/3) 



lim t -1 Tr ( n 

t— YOO 



e 2 



iA H + +iXfg dxh+(x) 



dx h+(x) 



At large times, the main contribution from the integral Jq dx h±(x) downstairs is that of points x 
lying far from and far from t. These points have approximately all the same contributions, and 
cover a domain of length diverging like t. Hence we can make the replacement f dxh+(x) i— > 
th+(t/2). We can then translate in order to obtain 







d{iX) 



/(A;/3) 



lim Tr ( n 

t— >co 



-W-iX/2)Hi 



H++iX f_^, 2 dx h+(x) 



Since we are evaluating the normalized average, with a complicated density matrix, of a local 

t /2 

operator h+(0), at large t, the operator J_ t , 2 dxh + (x) in the density matrix can be replaced by 
its limit h + (x) = H+. Hence, we find 



a 



8{iX) 



/(A;/3) = Tr n 



-W-iX)H^ 



h + (0) 



This means that 

a 



a(ix) 



F(X) = Tr n 



-(P t -iX)H+ 



h+(Q) 



Tr n 



-({3r+iX)H_ 



h-(0) 



Vfa-iXfir+iX 



(68) 



[^(0)]. 
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Integrating the result ([55]) . we obtain ([6T]) . 

We now prove the nontrivial steps above using our precise formalism. We start with three 
technical clarifications. 

First, in the precise definition of the operator (|65|) . we will require the notion of integral of 
local fields. The extension of the algebra V to include formal (algebraic) integral of local fields 
over finite intervals is completely straightforward, with in particular the supports being the 
(closed) integration regions. The action of the state un x n 2 on expressions containing integrals is 
immediately interpreted using Riemann integrals, which here always exist. 

Second, the operator (|65p will require an appropriate UV regularization. Given an e > 
small enough, we will define regularized states uj% and w| tat , acting on the linear space 
spanned by formal products of elements in V(M) <8> V^(M), by the requirement 

<stat [M x o) ■ ■ ■ Vn(x n )] = W ,stat [MVo) ' ' ' V n {y n )] \ yk=Xk+ikt (69) 

where on the right-hand side we make an analytic continuation (it is unique and well-defined). 
We will define <r| similarly. This requirement corresponds to the usual Euclidean time-ordering 
of QFT. As a result, ul,uj^ tat do not act on V ® V because the ordering is important, but the 
(weak) limits lim e ^.o+ w »i w Itat are zero on the ideal generated by the commutators, and, seen 
as acting on the quotients from the subspace where points are non-coincident in products, they 
give ^oj^stat respectively. We will see <J> , Ut, etc. as acting on V" ® in the natural way. 
Finally, we will use the large-distance "cluster property" (factorization property) of the state 

if 3 A £ Supp(C), B € Supp(£) : minflsc - y\ : x G A, y e B) oo. (70) 

This property is known from vertex operator algebras (it is a general property in QFT). 

We now provide the precise starting point for Proposition 16.11 The point + in (I57p corre- 
sponds to a limit process. Hence let us choose a a > small enough and define 

AE(t,a):= I dx (h + (-x) - h_(x)). (71) 

J a 

It is possible to show, tracing the steps below, that the large-t limit commutes with the a — > + 
limit; and that the large-t limit does not depend on a. Hence in the following we will make the 
simplifying assumption that |A| < 2a, in order to lighten the proof. Following the arguments in 
Section 14.51 the heuristic definition (|56p corresponds to 

Ad (e UE ) ^ Uf := $ D o (U\ /2 (8 U\ /2 ) o SJ 1 . (72) 

Note that the evolution operator Xjf preserves the initial state u a . Then, by a standard integral 
expression for the BCH formula, we may use the following expression for a regularization of O 
as defined in (|65|) . seen as an element of V* ® V': 

O i-> O a {\) := Pexp / dz{Wf (AE(t,a))) . (73) 
Jo 
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Here V denotes the ordering according to which integrands with values of z nearer to the upper 
limit A are positioned further to the right. Hence, our starting expressions are 

F vN (A) = lim hm t- 1 log(^[O a (A)]) 

F ind (A) = hm hint- 1 log K tat [O a (A)]). (74) 

We will see that the limits e — > + indeed exist. This fact indicates that the results are universal, 
and do not need renormalization. 

From (|74|) . we follow the lines of the heuristic proof presented above. The first step is to 
obtain 

F v nW = ^ind(A) = /(A; A) + /(-A; /?,) (75) 
where ^ 

/(-A;/3) := lim lim t~ l loga%\V exp [ dz ( -i [ dx Mx + z/2) ) 1 . (76) 

e->0+ <->oo Jo \ Ja J 

Consider F^. From the definition (|36p and from (|7ip . we see that for all t > a > we have 
$^ 1 (A£'(t, a)) = AE(t,a), that is, acts as the identity. Then, we find 

ft— z/2 ft-\-z/2 

(ul z/2 ®U^ z/2 )(AE{t,a)) = dxh+(-x)+ dxh.(x) (77) 

^ ' Ja-z/2 Ja+z/2 

and again $ G acts as the identity on this whenever \z\ < 2a. This implies that in (|73p we 
may replace XJf by U^_ z , 2 ® U^/ z ^ 2 . Thanks to (j4T)|) (and recalling that w s t a t = 0p u p r ), this also 
implies the first equation of ([ST]) (the first equality in ([75]) ). Using (j30|) . we immediately obtain 
the second term of the second equality of (|75p . For the first term, we obtain instead 



lim t 1 log cr^ PP exp / dz ( i dxt)(—x + z/2) 

J0 V J a 



We may use reflection invariance of the state ap, replacing —x + z/2 ^ x — z/2, and then make 
the change of integration variable z i— > —z in order to obtain /(A; /Jj). 

The second step is to take the A-derivative of (|76p . and make a translation in order to 
symmetrize the limits of the x-integrals (here we already put a = as it will be clear that this 
shift of the integral limits will not change the large-t result): 

df(X;P) . . ^[Pexpj- X dz {-iff /2 dx^{x + z/2)) f)(y - A/2)] 



— z — = lim lim t dy 

OX e^o+t^oc J_ t/2 a^[VexpJ- X dz(-iJ% 2 /2 dxi)(x + z/2y 

(78) 

By standard combinatoric arguments, the ratio in the y-integral is an infinite series of z- and 
x-integrals of connected averages of z-ordered factors: 



I'Sn'SS,'" 1 f_l dya >>t PeXP l ^(-if ' dxf,(x + z/lj) h(y- A/2)] coniiectrf . 



-t/2 JO x J-t/2 
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For convenience, we re-write the right-hand side of (|79p by extracting from the exponential 
the integration over x that runs from — oo to oo: 

,i/2 

lim lim t / dy ■ 

a% \ V e /o~ A dz d * "(-+-/2)) Jo* dz H-flT + f^ b(*+*/2)) hf _ A/2) i 

PL >\a I )\ connected 

Thanks to the cluster property ()70p . any connected average where fields are at positions that 
are far from each other tends to zero. Further, by standard results in CFT for states of the 
type ([23]) (or ([27]) ) such an average tends to zero exponentially fast for any /3 < oo. Hence: (i) 
taking the zeroth expansion term of the second exponential (i.e. replacing it by 1), the large-i 
limit of the second line (i.e. for fixed y) exists on every expansion term of the first exponential; 
and (ii) the large-t limit of the full integral over y exists on every expansion term of the first 
exponential and on every expansion term of order greater than zero of the second exponential. 
Thanks to the factor t , this implies that we may keep, in the large-t limit, only the zeroth 
expansion term of the second exponential. By translation invariance, we may fix y — A/2 to 0, 
and this gives 

ft/2 

j& & ^ L/2 dy °* [v eIt;X dz dxWx+z/2)) m] 

= \im + a^[Pe^ Xdz (- i ^ dxt,(x+z/2 ^ f)(0)] connected 

= M^[e^n^Wft(0)] comecte , (80) 

Finally we use the expression (|27|) of the state up from which the regularized state is 
defined, where we may replace dx t)(x + z/2) by its correspondent for finite R: 

lim lim a%. R \e iXlR K> xi,(x) {,(0)1 

e ^. +-R^oo P ' Hl n ^connected 

Here &p. R is the regularized version of the state (Jp-R which corresponds to (|26p . Thanks to the 
fact that the cluster limit is uniform in the large- R limit, the above replacement gives vanishing 
corrections as R — > oo. In cr^, R , we may use the standard algebraic construction of vertex 
operator algebras. Thanks to the e-regularization corresponding to the usual time-ordering, we 
see that the algebraic construction leads to the replacement 



i\ f R L 2 , n dx b(x) 



The normalization implied by n[-] results from the fact that we are evaluating connected averages 
in the initial expression. The result is now independent of e, because Lq is time-translation 
invariant. From (|27p we see that the above replacement corresponds to the shift j3 h- > (3 — iX. 
The result is then a one-point function, whose connected average is the usual average. Hence 
we find 

-i^%^ = ^- iX [m]- (8i) 

This is the equivalent of (|68p . and completes the proof. 
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7 Charge flow statistics and its large deviation function 



We now extend the results of previous sections to charge transfer. This happens when the 
conformal chiral algebra V contains an affine Kac-Moody sub-algebra on top of the Virasoro 
algebra @. One can then pick a u(l) affine sub-algebra of the Kac-Moody algebra and use it to 
measure what we call the charge. Generators j n , n € Z, of the u(l) affine algebra are chosen to 
be normalized according to [j n , j m ] = n<5 n+m; o. Let the charge number operator be the element 
91 := jo of V . Borrowing notations from previous sections, the it(l)-current is identified with 
K x ) = J-il( x ); normalized by the operator product expansion 

)( x )j(v) = 7 — — T2 1 + re S ular - ( 82 ) 
(x -yy 

Similarly as in previous sections, we assume having prepared the two halves of the system at 
equilibrium but with different temperatures and different chemical potentials (the later being 
coupled to the u(l) charge) and glued them at a very early time. We will construct the steady 
state of the total system, reached after a long waiting time, and compute the large deviation 
function of charge transfers. Proofs will be shorter than above since they are parallel to those 
used for energy transfer. 



7.1 Heuristics and the steady state 

The physical current is actually made of a pair of chiral currents: the left and right-moving 
currents j± which according to previous settings are identify as j+(x) := (g) 1 and j-(x) := 
1 <8> ](x). Their time evolutions are chiral, (dt ± d x )j± = 0, so that locally j±(x,t) = j{x =F t). 
The u(l) charge that we denote by iV to t is the integral of (j + +j—)(x) along the system interval: 

/•O- pR/2 

JVtot := / dx (j + + j_ ) (x) + / dx (j + + j_ ) (x) . 

J-R/2 J0+ 

It is conserved by the time evolution if j+(±R/2) = j_(±i?/2) and if j-)(x) is continuous at 
the contact point at x = 0. The dynamics before and after contact correspond to two different 
choices of boundary conditions: j + (0 ± ) = j-fi^) before contact and j±(0~) = j±(0 + ) after 
contact. 

We start with the uncoupled density matrix p = Pi® Pr & t time —t Q with 

p l>r : =tt[e-A.-^' r -^^' r )], 

where H l Q ' r are the energies and N l ,r the u(l) charge numbers on the left/right subsystems, 
naturally defined by restricting the integrals of the energy densities and the currents to the 
subsystem intervals, and N tot = N l Q + N'^. As usual, by duality this defines Gibbs states over 
the chiral algebras. 

9 More precisely, the chiral algebra V contains the semi-direct product of the Virasoro algebra by the affine 
Kac-Moody algebra since the Virasoro algebra acts naturally on the Kac-Moody algebra. 
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As above, the steady state /O s t a t is reached in the large time limit after the two subsystems 
have been coupled, and the 5-matrix intertwines p Q and p s tat: 

Pstat = S p Q S 1 . 

So we need to know the action of the S-matrix on the u(l)-currents. We have: 

S- 1 j + (x)S = j + (x), forx<0, (83) 

S^ 1 j+(x) S = j-(—x), for x > 0. 

The proof is exactly the same as for the chiral energy densities: from its definition, S := 
lim^^oo e~ toH e +t ° H ° in the large R limit, the 5-matrix acts on operators by first evolving them 
with the .ff-dynamics forward in time and then backward in time with the -ffo-dynamics, and 
these two only differs by the boundary condition at the contact point. Similarly, 

S~ 1 j-(x)S = j+(-x), forx<0, (84) 

j~(x) S = j-(x), for x > 0. 

Again, the left (right) movers end up in the right (left) part of the system, and according to 
the usual arguments, the steady state factorizes on left/right movers, 

where 0p lr ,n lr are the Gibbs states on the left/right chiral algebras at respective temperature 
Pi r and chemical potential /x; jr . 



7.2 Charge transfer 

The charge transfer that we will look at is that of half the difference of the u(l) charges on the 
left and right subsystems, Q := i(iV£ — N^). In the infinite volume Q is clearly infinite, but its 
time variation under the .ff-dynamics is finite. For i?>twe have 

AQ(t) := Q(t) -Q= f dxj + (x) - [ dxj_(x), 
J-t Jo+ 

which simply expresses the ballistic transport of charges. 

To get the large deviation function of the charge transfer we must evaluate the generating 
function Z t {v) := oj stat r e ^(Q+ A Q(*)) e -« i/ Qj ; which implicitly depends on f5^ r and pi yT . Following 
the by-now usual reasoning, Zt{y) factorizes into the product of two factors associated to both 
chiral sectors, and we may write 

Zt{v) = [0+] a^ r [O-] , with Q± = e^ N ^ iv So^j±(-) e ±l»N± 

Here N± refers to the charge number operators in the corresponding chiral sector (they arise 
because Qo act on left /right movers away from the origin as N + + iV_ ) , and we used translation 
invariance to replace the integral f® t by J *. Hence we have factorization of Zt, 

Z t iy) = zt(u; fa, m) zt(-u; fi r , /i r ) 
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for some single function Zt, which is defined purely in terms of the chiral algebra V , 

z t (u; ft, fj,) = [e-|^++-/o «faK-) e i** + ] . 

Assuming for a short while (but this will be proved in a few lines) that the following large- 
time limit exists, we define the charge transfer large deviation function F by 

F(v;Pi r ,ni r ) := lim -\ogZ t (v) 

t— >oo t 

Due to the factorization of Zt(u) we have: 

F(y;Pij.,nif.) = f(v;Pi,ni) + f(-v, P r , y r ), 

with f(v;ft,y) := lim^oo ~ \ogZt{v, ft, /i). By scaling invariance, f(u;ft,/j,) = ft~ x h{v;ft\x). We 
now claim 

Proposition 7.1 flj T/ie derivative of the chiral large deviation function f is related to the 
current one-point function in the Gibbs state but at a shifted complex chemical potential. Ex- 
plicitly 

d 

- ^/(^ A A 1 ) = <T l3,ii+iuf}-i [j(0)] , (85) 

(2) As a consequence, f{y; ft, fx) = [(fty + iv) 2 — (ft y) 2 ]/ (Sir ft), for unitary theory. 

(3) The charge transfer large deviation function is 

v(ll .a „ \ (Pm + iv) 2 - [ftm? , {ftrjh- - W? ~ {ftr^r? ^ 
* (y, fti r , M r) = — - 1 — — , (8b J 

87rp r 



and the statistics of charge transfer is gaussian with mean ^'^^ and covariance ^' 4 ^ r - . 

The proof of this proposition is based on similar arguments as that for energy transfer plus an 
input from conformal field theory techniques which allows us to evaluate the current one-point 
function. 

Proof. 

First we take the derivative of log zt(v; ft, y) w.r.t. iv. This has two effects: it inserts J dxj(x) 
inside the expectation value and it divides this expectation value by Ztiy; ft, By translation 
invariance, we may replace the large t limit of the integral by the insertion of t)(t/2). The last 
factor t combines with the factor 1/t in the definition of /, hence 

a^[e-t^++-/o^)W j (t/2) e +i^ + ] 
—ia v j[v; ft, a) = hm — — ; - t : -. 

'^'^ j- e -^JV++w;*d IB K*) e +^+] 
Before taking the large time limit we manipulate this expression: using translation invariance 

t /2 

we centred the inserted current at instead at t/2, and we then convert the integral /_ t / 2 dx)(x) 
into N + — ( f^ 2 + f_^ 2 )dyj(y). Recalling the definition of the Gibbs state 073^, this yields to 

-id v f{v;P,n) = lim L . H )\ _ 



t— ¥00 
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Now, as for energy transfer, the perturbative expansion of this ratio in a power series of v 
involves connected correlation functions of products of the current operator j(0) times series of 
multiple ordered integrals of j(y). Since the currents are located at distance at least t/2 from 
the origin, the clustering property of connected correlation functions ensures that we can safely 
take the large t limit. All contributions go to zero except the zero order in the perturbative 
expansion, and this amounts to take directly the large t limit in the above equation. Thus 

-id u f{v; (3, y) = <Jp^ +iv p-i [j(0)] . 

Actually, the proof can be slightly simplified using the fact that N + commutes with the current 
operator ){y) in the pure u(l) case that we are presently dealing with, but we need to keep this 
more general argument for next section. 

We proved point (1) of the proposition. To prove point (2) we need to compute <jp t(Jt [j(0)] 
and then shifts /3/i into /3fi + iv. This is done using conformal field theory techniques and 
modular transformation. Let x[r, z] be affine u(l) character: xl 7 "?- 2 ] = Tr(e 27 " T '' io_c / 24 * > e 47r * 2: ' ) 
with normalization j(z)j(0) ~ l/(z 2 ). Since according to Section EOl H + = ^-Lq at finite size, 
we need to set r = i/3/R and z = — i/3///(47r), so that 

°pM°)) = J™o^?|l l0gx[ * /3/ ^'"^ /(47r)] ' 

by translation invariance. To evaluate the large size limit we implement a modular transforma- 
tion. Modular property of characters of affine algebras are well known. They form representa- 
tions of the modular group, and we may write 

XM = e~ l2 ^ lT "£S pXp [-1/t,-z/t}. 
p 

where the sum is over a set of affine u(l) modules, including the identity module V, and S p are 
coefficient of the so-called modular S- matrix [TU]. For unitary theory, the leading contribution 
for R large comes from the identity module in the sum. This gives 

WXO)] =/V(4?r). 
Integrating —id v f = a a u+ivp- 1 D(0)j yields the formula for /. 

8 Energy and charge flow statistics and its large deviation func- 
tion 

We now look at the large deviation function combining energy and charge transfers associated 
to the two commuting observables E := \{H r — H l a ) and Q := \{N^ — N l Q ). As above we have 
to compute the generating function Z t {\,v) := w stat [e^W e lXE ^ eT %v Q e~ lXE ], and the large 
deviation function 

F(\,v;/3i r ,(n r ) := lim - log Z t (X, v). 
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As before, Zt factorizes in left /right expectations so that the large deviation function reads 

F(\,v; fl,r,w, r ) = f(\,v,/3 h fj,i) + /(-A, -v;f3 r ,(j, r ). (87) 

The function f(\,v; ft, n) is the chiral large deviation function which is computed using data 
from the chiral algebra V only: 

/(A, v\ (3, n) = hm \ log [0\ e^ N + ■ 0\ e^ XH +] 

t—t-oo t 

with op^ the Gibbs state over V at temperature 1//5 and chemical potential //, and 

0\ = e -i XH ++ iX fo dxt ^( x \ 0\ = e ~i uN++iv fo dx '^ x \ 

Let us now compute /. We follow the same strategy as in previous section, and we shall be 
more brief. As before, we take derivatives w.r.t. A and v. Using translation invariance to replace 
the integrals of the energy or charge densities by the insertion of the corresponding operator at 
the middle point t/2, we get 

.« tn o s r ^A°\^ vN+ -oi^m^] 

-iOxf{X,v,p,H) = hm 



— lO u j{\,V\p,Ll) = lim ; ; . 

This ratio of expectations only involve connected correlations functions. So, by the usual argu- 
ment this garanties that we can take the naive large time limit which amounts to replace O l e by 
e 2 XH + and O l by e? vN+ . As a consequence, we get 

Proposition 8.1 (1) We have: 

— id\f(\, v;p~,fj,) = a ? iX ^[W)), (88) 

-id u f(\,u;/3,[i) = <Tp_ iK §0v [j(0)]. (89) 

(2) For unitary theory, we have 

/(A, u; p, n) = - (— --) + ^ _ , A) - — . (90) 

The large deviation function for energy and charge transfers is given by: F(\,v; $ r , /i/ >r ) = 
/(A, i/; A, Mi) + /(-A, A-,Mr)- 

Proof. 

We already argued for point (1). To get point (2) we have to evaluate both one point func- 
tions (|88l89p . This is done as previously, using modular property of characters and dp log x = 
-RapAW) ~ and d^logx = PR crp^ [j (0)] . For unitary theory, we get: 
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Notice the extra term in <7g )/U [f)(0)] for non zero chemical potential \i. This terms ensures that 
the two differential equations (|88|89|) are compatible. Integrating them with boundary condition 
/(0, 0; P,n) = yields formula (l90j) . 

Formula (|9U|) indicates that charge and energy flows are not statistically independent, even 
at large time, and this is a direct consequence of eq. (|9ip for the one-point chiral energy density at 
nonzero chemical potential. As a consequence, at identical temperatures but different chemical 
potentials, the mean energy current does not vanish (Je) A^=o,A^o = ^(a 4 ? ~ A*r)- 

9 Conclusion and discussion 

We have constructed non-equilibrium steady states in conformal field theories and derived the 
associated large deviation functions. It is worth noticing the universal character of these states. 
As a consequence the large deviation functions for energy and charge transfers are very univer- 
sal, depending only on universal constant (h, ks) and on the temperatures T^ r and chemical 
potentials /i;. r . Although we used a conformal field theory approach, we provided arguments 
in favour of the universality of our construction in the low energy regime. Remark also that 
we would have got identical steady states if we would have prepared the two infinite left /right 
subsystems at different temperatures and connected them with some finite-length subsystem 
initially in another distribution. It should also be noticed that naive Fourier's law for the energy 
current as a function of the temperature gradient is broken in our model 0. Indeed, in our 
model, the temperature profile, say specified by the local energy density, is flat, although the 
energy current is non zero (this may be related to the zooming prescription that we alluded 
to in Section [2]) . As we shall explain [5] , the key equations (|88l89p relating derivatives of the 
large deviation functions to the linear response theory but at shifted temperatures and chemi- 
cal potentials fit into consequences of PT-symmetry, as does the reflection-less Lesovik-Levitov 
formula for charge transfer in mesoscopic free electronic systems. Generalisations of our con- 
struction to include non-topological defects, or massive integrable perturbations, or to higher 
dimensions should provide valuable information on non-equilibrium quantum physics. 

Acknowledgements: We thank J. Cardy and M. J. Bhaseen for interesting discussions. 
This work was in part supported by ANR contract ANR-2010-BLANC-0414. 

A Two-time energy measurement in the steady state 

Here we briefly present the derivation of the large deviation function according to the measure- 
ment scheme described by eqs. (|63!64p . According to (f64"|) . take 

O' = e -<!-«)£ e iM£+A£(t)) e -i(!+«)E_ (92) 

10 In low dimension one naturally expects a breakdown of the Fourier's law with energy current scaling as L~ a 
for some exponent a if the two reservoirs are at a distance L apart. 
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We follow the lines of reasoning of Section \6. 11 Denning AE(t, a) as in (|7ip and using the BCH 
formula, we may represent the regularization of O' by 



O' ^ O a (X,u) := U*_x/2 (Vexp j\z (iU? (AE(t,a))) 



(93) 



with Uf defined in (|72j) . In order to describe the two-time measurement, recall that we need 



to integrate over the variable u, see ()63j) . In the infinite- length limit, the energy spectrum 
is continuous, hence the integration range is (— r, r) with r — > oo. The precise proportionality 
factor in (|63|) is not important, except for the fact that it becomes proportional to 1/r as r —¥ oo. 
Hence, our starting expressions are 

F(X) = Km i^log (j^r' 1 J duu> sta t[O e , a (\,u)]J (94) 

Factorisation of w sta t implies 

F(A) = /(A;A) + /(-A;/3,) 



where 



f(-X;(3) := lira t _1 log<r^[ ^ ex P^ ^ fj(a; + z/2)^ 



(95) 



The main point is to notice that the integrand inside the u integral in the first equation of (|74p 
does not vanish as \u\ —¥ oo, but rather converges to a constant. Hence, the A — > oo limit is 
given by this constant. In order to evaluate this constant, notice that 

U*- X/2+z (^E(t,a)) (96) 



We then use 



x — — + — — - J— + — — — + u large and positive 

-x + ^- ^ + ^J-/i + f-x-^ + ^-^jj u large and negative. 



O a (X,u) := (vexpj dz (i U*_ x/2+z (AE(t, a)))) 



and for large |u|, we observe that both cases u > and it < in (|96p involve two terms largely 
separated, so that we can make use of the factorization property (fTOj) of Wgtat = &/3 lt i3 r - Using 
translation invariance the u-dependence of the result disappears. We can also use translation 
invariance to take away the explicit A-dependence in the arguments of h±, and then we may use 
the change of variable x ^ t -\- a — x followed by a translation by ±(t + a) in order to change 
the sign of x. This allows to bring each factor to a form involving the path-ordered exponential 
in (75]). We then find 



limr" 1 duu st!lt \O ea (X,u)] ^ e *(/(^)+/(-A;A-)) + e *(/(A;A)+/(-A ; A)) as t _> ^ (97) 
This shows that the differences between -F(A) and -F(A). 
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